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Abstract 



A measurement of the absolute fluorescence yield of the 337 nm nitrogen band, 
relevant to ultra-high energy cosmic ray (UHECR) detectors, is reported. Two in- 
dependent calibrations of the fluorescence emission induced by a 120 GeV proton 
beam were employed: Cherenkov light from the beam particle and calibrated light 
from a nitrogen laser. The fluorescence yield in air at a pressure of 1013 hPa and 
temperature of 293 K was found to be Y337 = 5.61±0.06 S £ a t ±0.21 syst photons/MeV. 
When compared to the fluorescence yield currently used by UHECR experiments, 
this measurement improves the uncertainty by a factor of three, and has a significant 
impact on the determination of the energy scale of the cosmic ray spectrum. 

Key words: Nitrogen Fluorescence Yield, Air Fluorescence Detection, Ultra-High 
Energy Cosmic Rays 

PACS: , 96.50.S-, 96.50.sb, 96.50.sd, 32.50.+d, 33.50.-j, 34.50.Fa, 34.50.Gb 



1 Introduction 



A well established technique for the detection of Ultra-High Energy (> 10 18 eV) 
Cosmic Rays (UHECRs) - first successfully employed by the Fly's Eye [1] and 
HiRes [2] experiments - is based on nitrogen fluorescence light emission induced 
by Extensive Air Showers (EAS). Excitation of atmospheric nitrogen by EAS 
charged particles results in fluorescence emission, mostly in the wavelength 
range between 300 and 430 nm. This UV light is measured as a function of 
time and incoming direction by photomultiplier cameras at the focus of large 
(few m 2 ) mirrors. Fluorescence telescopes measure the longitudinal EAS devel- 
opment in the atmosphere, which provides unique information on the primary 
cosmic ray's type and a calorimetric measurement of its energy. 

The fluorescence light yield along the EAS depends on the air pressure, tem- 
perature and humidity at the point of emission. In addition, wavelength- 
dependent atmospheric attenuation affects the light intensity reaching the tele- 
scope. Thus, the intensities of the fluorescence bands must be known for atmo- 
spheric conditions corresponding to the EAS development in the atmosphere, 
which ranges between about 2 km and 15 km above sea level. Early measure- 
ments of the fluorescence yield include those with low-energy stopped-particles 
in air by Bunner [3] and with electrons in air by Davidson and O'Neil [I]. A 
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number of other experiments have made measurements of the fluorescence 
light yield pertinent to EAS [6] [7] [8] [9] [10]. 

The AIRFLY (AIR FLuorescence Yield) Collaboration has carried out an ex- 
tensive program of measurements to significantly improve the precision on the 
fluorescence light yield. The fluorescence emission was studied as a function 
of the kinetic energy, ranging from keV to GeV, of particle beams at several 
accelerators [11]. The relative intensities of 34 fluorescence bands in the wave- 
length range from 284 to 429 nm, together with their pressure dependence, 
were reported in [12]. The temperature and humidity dependence of the main 
fluorescence bands was also measured [13] . These measurements have provided 
the most complete and consistent set of fluorescence yield data for UHECR 
calibration, establishing the relative spectral line dependence on atmospheric 
parameters. In this paper, a precise measurement of the absolute fluorescence 
yield - namely the number of photons emitted per MeV of energy deposited - 
is reported. This result will impact the current generation of UHECRs experi- 
ments, the Pierre Auger Observatory [14] and the Telescope Array |15j . as well 
as a future large observatory or the proposed space-based JEM-EUSO [16], 
which rely on fluorescence detection for their absolute energy scale. In fact, 
the fluorescence yield is presently a major contribution to the systematic un- 
certainty of the cosmic ray energy (14% for the Pierre Auger Observatory [T7] 
and 11% for the Telescope Array [18], with both experiments quoting a total 
systematic uncertainty of 22%). 

In a standard approach to a photon yield measurement, the absolute quantum 
efficiency of the photon detector ultimately limits the systematic uncertainty. 
AIRFLY follows a different approach, first presented in [19], based on the 
comparison of the fluorescence yield to Cherenkov light emitted in the same 
experimental setup. Measurements are performed with a narrow band optical 
filter centered around the main nitrogen fluorescence line at 337 nm. Since 
the fluorescence yield is measured relative to the known Cherenkov yield, 
the absolute value of the photon detection efficiency is not essential, with a 
considerable reduction of the systematic uncertainties. In this work, several 
improvements have been implemented with respect to the preliminary test 
in [19], most notably the use of a single particle beam at the Test Beam Facil- 
ity of the Fermi National Accelerator Laboratory (FNAL) and of an integrat- 
ing sphere for optimal light collection. Measurements were mostly performed 
with 120 GeV proton© In addition, a measurement of the fluorescence yield 



Notice that, at these energies, the nature of the primary particle does not affect 
the fluorescence emission of the 337 nm band, which is induced by the secondary 
electrons produced in the gas |20| . In fact, the cross section for electron impact 
excitation of the C 3 II U state of the nitrogen molecule (responsible for the 337 nm 
emission) becomes completely negligible for energies above 10 keV, and direct exci- 
tation by protons is also negligible [2Tj . 
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was performed which employed a nitrogen laser as a calibrated light source. 
When the two methods are combined, a total uncertainty of 4% is achieved. 
A preliminary account of this measurement was presented in |22j. 




This paper is organized as follows. Section [2] presents the apparatus and ex- 
perimental methods. Simulations are described in Section [3j Data analysis 
is presented in Section HI The procedure for the fluorescence yield measure- 
ment is described in Section [5j Results of the fluorescence yield measurement, 
including systematic uncertainties, with the Cherenkov calibration and with 
the laser calibration are given in Sections [6] and [7J respectively. The combined 
measurement is presented in Section [HI and conclusions are drawn in Section [HI 



2 Apparatus and experimental methods 

In this Section, the AIRFLY apparatus is presented. Details of the integrating 
sphere, a fundamental component of the experiment, are given. The proce- 
dures for the measurement of fluorescence and Cherenkov light are illustrated. 
The laser setup, which provided a calibration light source for the fluorescence 
measurement, is also described. Lastly, the measurement of the 337 nm fluo- 
rescence line yield in nitrogen to that in air is presented. 

2.1 The AIRFLY apparatus at the FNAL Test Beam Facility 



The measurements were performed at the FNAL Test Beam Facility! 2 1, which 
provided a primary 120 GeV proton beam from the Main Injector and sec- 
ondary beams of muons, pions, electrons and positrons of lower energy. The 
proton beam was selected in the majority of the measurements for its larger 
multiplicity, smaller beam profile and negligible contamination. Every minute, 
the beam is spilled in the experimental area during 4 s, with particles grouped 
in trains of bunches. Trains are separated in time by 10 yus, and bunches within 
each train are separated by 19 ns. Typical conditions for data taking with the 
proton beam were 30 bunches per train with a multiplicity of 2 • 10 5 particles 
per spill. The beam profile, monitored by wire chambers placed before and af- 
ter the AIRFLY apparatus, was typically 3 mm x 4 mm (r.m.s.). The layout of 
the apparatus is shown in Fig. [TJ A scintillator disk counter (10 mm diameter 
and 3 mm thickness) tagged the beam at the entrance of the AIRFLY cham- 
ber. A Cherenkov counter - a UV-transparent acrylic cylindrical rod of 10 mm 
diameter and 30 mm long - tagged the beam at the fluorescence chamber's 
exit. Cherenkov light produced in this counter by the beam particles provided 

2 http://www-ppd.fnal.gov/FTBF/ 
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Fig. 1. The AIRFLY experimental setup: 1) stainless steel chamber; 2) integrating 
sphere; 3) top port; 4) plug of exit port; 5) 337 nm interference filter; 6) acceptance 
apertures; 7) PMT; 8) extension for beam particle path; 9) scintillator disk counter; 
10) Cherenkov counter, with acrylic rod and photomultiplier; 11) veto counters. 




a signal of few ns time resolution and very good single particle resolution. Two 
large 150 mm x 150 mm scintillator counters, with a central hole of 10 mm 
diameter for the beam, vetoed off-axis particles. The fluorescence chamber was 
placed between the two veto pads. The gas chamber was a 3 mm thick stainless 
steel cube of 200 mm side length. Flanges were provided for beam windows, 
shutters, gauges, gas inlet and vacuum. The beam entrance window, made of 
0.1 mm thick aluminum, was installed on an upstream extension of the cham- 
ber to provide an additional length of 18 cm for Cherenkov light production. 
The exit window, also of 0.1 mm thick aluminum, was mounted on an exten- 
sion to the downstream side of the chamber. An integrating sphere located at 
the center of the gas chamber (Section 12.21) collected the light produced by 
the beam in the chamber. Light out of the integrated sphere passed through 
a 337 nm interference filter (bk Interferenzoptik, model bk-337 1-10-C2, 70% 
peak transmission) that served at the same time as a chamber window. A 
mechanical optical shutter was incorporated with the filter. The shutter could 
be closed and opened remotely to alternate background measurements during 
data taking. A photomultiplier tube (Hamamatsu, model H7195P), hereafter 
called PMT, with good single photoelectron resolution was used for photon 
detection. The optical field of view was defined by a 40 mm diameter accep- 
tance cylinder placed between the integrating sphere's port and the filter, and 
by circular apertures of the same size placed in front of the PMT photocath- 



5 



ode. The chamber was internally covered with black UV-absorbing material to 
avoid stray light. An additional scintillator counter vetoing off- axis particles 
that could hit the PMT and a second H7195P photomultiplier with a closed 
shutter for background monitoring, both placed on the side of the PMT, are 
not shown in Fig.[TJ The chamber and the counters were mounted on an optical 
breadboard for precise alignment. The breadboard was placed on a remotely 
controlled position table, which allowed the alignment of the apparatus with 
the beam by maximizing the rates of the entrance scintillator disk and exit- 
ing Cherenkov counters. A laser system (Section 12. 3p was routinely inserted 
in the apparatus for calibration purposes. A remotely controlled gas handling 
and vacuum system was installed. The chamber pressure, temperature and 
humidity were measured. Typical pressure and temperature ranges during the 
fluorescence measurements were 1013 ± 10 hPa and 293 ± 2 K. The same high 
purity dry gases, 99.9999% pure nitrogen and an air-like mixture (79% ni- 
trogen, 21% oxygen), of our previous measurements [T2] P2] were used. Pure 
helium was used for background measurements. 

The data acquisition was based on a VME system. An event trigger was formed 
from the coincidence of a train trigger gate, issued in time with the arrival of 
each train of bunches, and a single particle trigger gate from a beam moni- 
toring scintillator counter. The Test Beam Facility provided both signals. A 
single detected particle within a train caused the signals from the apparatus 
scintillator counters and PMTs to be digitized by a 12-bit 500 MHz Flash 
ADC (model SIS3350, Struck Innovative Systeme). For each event trigger, 
600 samples equivalent to 1.2 /is - thus containing the entire train of bunches 
- were saved in the FADC memory. The train trigger gate was used as time 
reference for the trigger signal, and was found to be stable at the level of 
1 ns. As an example, the ADC trace of the Cherenkov counter for one event 
trigger is shown in Fig. 12(a), presenting two fast pulses produced by particles 
in two separate bunches. The ADC trace averaged over all event triggers in 
one spill is shown in Fig. [2(b), with the time structure of the bunches clearly 
distinguishable. 

The module's internal memory of 128 MSamples/channel was used to store 
the data during the 4 s spill without dead time. Data for up to 1 GByte 
(the exact number depending on the number of triggers during the spill) were 
read out and saved to disk in about 40 s between spills. While demanding 
in terms of data transfer and storage, this data acquisition strategy provided 
accurate information for each triggered train, allowing flexibility, redundancy 
and detailed systematic studies in the subsequent offline analysis. Data tak- 
ing was organized in runs integrating 30 spills for about 30 min duration, to 
guarantee stable conditions. A fast analysis of a sizable fraction of the data 
was performed online, providing several monitoring variables like rates, time 
distributions and ADC integral distributions. A total of 7 TByte of data was 
collected during two dedicated test beam periods. 
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Fig. 2. Digitized ADC trace of the Cherenkov counter: a) trace recorded for one 
event trigger, with two fast pulses; b) ADC trace averaged over all event triggers in 
one spill, showing the bunch structure of the beam. 



2.2 Integrating sphere and data taking configurations 



The integrating sphere - a hollow sphere with an internal surface highly diffu- 
sive in the UV - was a fundamental component of the AIRFLY apparatus. The 
sphere collected the light over almost the 47r solid angle. Since only about 0.2 
fluorescence photons of the 337 nm band were produced from a single beam 
particle in the chamber, light collection over a large solid angle maximized the 
signal over background and guaranteed a reasonable data taking time. Also, 
the light is made isotropic and unpolarized after several diffusion bounces on 
the sphere's surface, producing at the detection port a Lambertian light out- 
put which is practically independent of the original light distribution. Thus, 
the isotropic fluorescence light and the highly directional Cherenkov or laser 
light present an identical distribution at the entrance of the 337 nm filter, 
greatly reducing systematic uncertainties. 

The integrating sphere was formed by two hollow aluminum half-spheres of 
15.2 cm diameter and 3 mm thickness. Four ports of 38.1 mm diameter were 
machined on the sphere: two of them along the beam direction (beam entrance 
and exit), one at the side, orthogonal to the beam direction (photon detection), 
and the last port at the sphere's top. The sphere's internal surface was covered 
with a 0.5 mm thick GORE™ DRP® 123 reflector. This proprietary material, 
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made of a micro porous structure of PTFE, has a very high diffusive reflectance 
(~ 99.5%) around 337 nm, providing a sphere throughput ~ 40% higher than 
that of a commercial sphere made of granular PTFE. The thin GORE™ 
DRP® reflector, when compared to the many mm thick PTFE walls of a 
commercial sphere, minimized Cherenkov light produced by particles passing 
through the diffusive material. Two remotely operated shutters could close 
the downstream exit and the top ports of the sphere, while the entrance and 
the side ports were always open. The port plugs on the shutters were made of 
the same GORE™ DRP® reflector, hereafter called diffuser for the sake of 
simplicity. 

Two different configurations of the chamber were used to detect either fluo- 
rescence or Cherenkov light. In the Fluorescence configuration, the exit port 
of the integrating sphere was open, while the top port was closed. In this way, 
the highly directional Cherenkov light, emitted along the particle path in the 
chamber, does not hit the sphere's diffusive surface. After exiting the sphere 
through the downstream port, the light is absorbed by the black inner coating 
of the chamber, acting as a Cherenkov light dump. Thus, only fluorescence 
light emitted along the particle path is collected by the integrating sphere 
and measured by the PMT viewing the sphere's side port. In the Cherenkov 
configuration, the top port was open, while the exit port was closed. With this 
arrangement, the Cherenkov light is diffused back into the sphere by the exit 
port plug, the integrating sphere collecting then both the Cherenkov and the 
fluorescence light. Notice that, by opening or closing the top port, the total 
reflecting surface of the integrating sphere is kept the same in both configura- 
tions - always three open ports - allowing for an accurate comparison of the 
fluorescence and Cherenkov measurements. 

Measurements were also performed with aluminized mylar mirrors replacing 
the diffuser plugs at the top and exit ports. The mirror reflective surface was 
placed at an angle of 160° with respect to the beam direction, thus reflecting 
the Cherenkov light back into the sphere. The use of a thin mirror rather than a 
diffuser plug eliminated a major background - Cherenkov photons produced by 
beam particles in the PTFE of the exit port plug. The top port was also fitted 
with a mirror, to keep the sphere port's symmetry between the Fluorescence 
and Cherenkov configurations. 



2.3 UV Laser calibration 

The laser calibration system consisted of a 337 nm nitrogen laser (Spectra- 
Physics, model VSL-337ND-S), producing light pulses of ~ 3 ns width at up 
to 20 Hz repetition rate, and a PTFE integrating sphere (SphereOptics, model 
SPH-6Z-4) referred to as laser calibration sphere in the following. A sketch of 
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Fig. 3. The laser calibration system integrated in the AIRFLY setup: 8) and 10) 
collimators; 9) laser calibration sphere; 11) optical attenuators; 12) 337 nm nitrogen 
laser; 13) PMTi aser , to monitor laser intensity. For all other components, see caption 
of Fig. CD 

the laser system during a calibration run at the Test Beam Facility is shown in 
Fig. [3j The laser light was attenuated by a set of filters and collimators before 
entering the laser calibration sphere. A photomultiplier tube, PMTi aser , which 
detected light reflected by an attenuation filter at « 180° from the laser beam, 
was used to provide a fast trigger signal of the laser pulse. Light exited the 
laser calibration sphere through a 6 mm diameter hole at the end of a 15.5 cm 
long collimator tube mounted on the sphere's exit port. The efficiency of the 
laser calibration sphere was determined from the ratio of the light pulse energy 
measured at the exit of the collimator tube to that measured at the entrance 
of the sphere. A silicon energy probe (LaserProbe inc., model RjP-465) with 
fj sensitivity was used to measure the laser pulse energy. The efficiency of 
the laser calibration sphere was measured to be (2.49 ± 0.02) • 10~ 6 , where 
the 0.8% uncertainty takes into account the probe linearity (0.5%) and the 
reproducibility of the result (0.6%). In a laser calibration run, the system 
was mounted at the front of the AIRFLY chamber, with the collimator tube 
facing the entrance port of the integrating sphere (Fig. [3]). The laser pulse 
energy at the entrance of the laser calibration sphere was ~ 5 pJ, measured 
with 5% uncertainty by the RjP-465 probe, which gave a calibrated beam of 
~ 20 photons per pulse entering the AIRFLY sphere. For each laser pulse, the 
PMT signal was digitized by the same electronics chain of the fluorescence 
measurements (Section 12. ip and a single photoelectron signal was searched for 
(Section Hj). A calibration constant, S^ser, was then obtained by counting the 
number of detected photoelectrons per 337 nm photons entering the AIRFLY 
sphere (shortened as pp in the following). A typical measurement yielded 
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'S'laser ~ H ■ 10~ 4 pp. The laser beam energy was measured by the silicon 
probe at the beginning and at the end of each 15 min calibration run. The 
stability of the laser pulse energy during the calibration measurement was 
monitored by recording the signal from PMT laser , and found to be typically 
better than 1%. Several laser calibration runs were taken each day at the end 
of the fluorescence measurements. 



2.4 Nitrogen to air fluorescence ratio 



With the beam multiplicity available at the Test Beam Facility, the detected 
rate of fluorescence photons in air was very small. For a precise measurement, 
long data taking runs would have been necessary, which were unpractical in 
the context of a test beam. Thus, most of the measurements were performed 
with pure nitrogen, where the fluorescence yield in the 337 nm band is almost 
ten times larger. The air fluorescence yield was then obtained by using the 
independently measured ratio of the 337 nm fluorescence in pure nitrogen to 
that in air, r^ 2 . 

AIRFLY measured the nitrogen to air fluorescence ratio as a function of pres- 
sure in [12], obtaining tn 2 = 7.38±0.08 at 1000 hPa. The ratio could be slightly 
different in the present setup, which included the integrating sphere and a new 
337 nm interference filter. Thus, a dedicated measurement was performed with 
a 18 /iCi 241 Am source. The radioactive source was placed at 3 cm distance 
from the entrance port of the sphere. Alpha particles from 241 Am have a range 
of about 4 cm in air, thus depositing most of their energy close to the entrance 
port. The air fluorescence spectrum induced by alpha particle excitation [23] is 
consistent with the AIRFLY measurement with electrons [T2]. Measurements 
of the PMT single photoelectron rates were then performed in pure nitrogen 
and air. In order to keep the same path length of the alphas in the gas, and 
thus the same spatial distribution of the fluorescence emission, measurements 
in nitrogen and air were taken at slightly different pressures to compensate 
for the gas densities. From these measurements, a ratio 

r N2 = 7.45 ± 0.07 (1) 



at 1000 hPa is derived, in very good agreement with our previous measure- 
ment [12]. 
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3 Monte Carlo simulations 



A full simulation of the AIRFLY apparatus was performed in the Geant4 [25] 
framework. In this Section, details of the photon detection components - the 
integrating sphere, the 337 nm filter and the PMT - are given. Also, the 
implementation of the fluorescence and Cherenkov yields in the simulation 
is described, and studies of the energy deposited by protons in gas are pre- 
sented. Finally, results of the simulation relevant for the fluorescence yield 
measurement are given. 



3.1 Photon detection efficiency 

The diffuser coating the inner surface and the port plugs of the integrating 
sphere was simulated as a 0.5 mm thick PTFE layer, diffusing incident light 
with reflectivity e^g and a Lambertian angular distribution. The value of e^s 
was chosen to reproduce the sphere's throughput measured in a dedicated ex- 
periment. A 337 nm laser beam was shot into the sphere, and the laser light 
power was measured at the sphere's entrance and at 11.5 cm distance from 
the side port. The ratio of these two measurements provided a benchmark 
efficiency of the sphere to be reproduced by the simulation. The efficiency, 
6.04 • 10~ 4 , measured with a 1% uncertainty which includes the photodiode 
non-linearity and the reproducibility of the result, was reproduced in a full 
simulation of this experiment for e^m = 99.43%. The sphere with aluminized 
mirror plugs was characterized by a similar set of measurements, which were 
reproduced by introducing in the simulation a value of 87.2% for the reflec- 
tivity of the aluminized mirror at 337 nm. 

The wavelength dependence of the integrating sphere's throughput was also 
included in the simulation. For this purpose, the reflectance of the diffuser 
was measured as a function of wavelength with a standard reflectance setup 
(integrating sphere AvaSphere-50 and spectrometer AvaSpec-2048). A weight 
was then implemented in the simulation for each diffusion process occurring 
inside the sphere. The aluminized mylar mirror reflectivity was also measured, 
and found to be constant within the bandwidth of the 337 nm interference 
filter. Since photons bounce on average about 20 times on the diffuser inside 
the sphere before exiting from the side port, the wavelength dependence of 
the sphere's efficiency is determined by the diffuser properties, and it did not 
change significantly with the mirror plugs. 

The transmittance of the 337 nm interference filter was measured as a func- 
tion of wavelength and angle of incidence of the light with a Lambda EZ210 
UV/Vis spectrophotometer. This dependence must be included in the simula- 
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tion, given the broad spectrum of Cherenkov light and the strong dependence 
of interference filter's transmission on angle of incidence. Measurements in 
steps of 0.1 nm and for 13 incident angles between 0° and 20° (the maxi- 
mum angle of incidence in the AIRFLY geometry) were performed. Data were 
linearly interpolated to cover all wavelengths and angles in the simulation. 

The quantum efficiency of the PMT was provided by the manufacturer as 
a function of wavelength. The absolute detection efficiency of the PMT was 
implemented in the simulation as the product of the quantum efficiency by 
a collection efficiency of 85% (a typical value for this type of photomultipli- 
ers according to the manufacturer). The absolute efficiency of the PMT as 
a function of wavelength was also independently measured with a dedicated 
setup, to be reported elsewhere [20], yielding a result consistent with the 10% 
systematic uncertainty quoted by the manufacturer. 



3.2 Photon yields and energy deposit 

In the Geant4 simulation, fluorescence photons were produced proportional 
to the energy deposited in a given air volume, and then emitted isotropically. 
The proportionality of the fluorescence yield to the energy deposit has been 
experimentally verified [7] [8] [11] and is justified by the physics processes 
occurring at the microscopic level [20]. Photons were sampled from the air 
fluorescence spectrum measured by AIRFLY [12] , with a nominal yield for the 
337 nm band of 5.15 photons/MeV of energy deposit in air at a pressure of 
1013 hPa and temperature of 293 K. 

Particles, both primaries and secondaries, are individually tracked in the sim- 
ulation and energy is deposited along their path according to standard electro- 
magnetic and hadronic processes. To speed up simulations, when the particle's 
range becomes shorter than one mm, tracking is stopped and the particle's 
energy is deposited locally. This procedure provides a realistic simulation of 
the spatial distribution of the energy deposit (and consequently of the fluores- 
cence emission), which is very important to avoid biases in the measurement 
of the fluorescence yield [2"T] . 

Dedicated simulations were performed to check the Geant4 model for the 
energy deposit. Results are shown in Fig. HI Reference calculations of the 
proton energy loss in nitrogen from the Particle Data Group [2B] and NIST [21] 
are accurately reproduced by the simulation. Also, simulations of the energy 
deposit in the integrating sphere's nitrogen volume are well matched by the 
restricted proton energy loss [28] calculated for a maximum energy transfer 
of 100 keV to secondary particles. This result is consistent with the path of 
S rays being limited by the size of the sphere, since electrons of 100 keV kinetic 
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Fig. 4. Simulation of proton energy deposit in nitrogen. Standard energy loss param- 
eterizations from the Particle Data Group and NIST are compared with that im- 
plemented in Geant4 (full lines, differences undistinguishable on this scale). Closed 
dots represent the energy deposited by a proton in the integrating sphere's volume 
according to Geant4; superimposed as a dashed line is the restricted proton energy 
loss calculated for a maximum energy transfer of 100 keV. 

energy have a range of 14 cm in nitrogen. Notice that a 120 GeV proton has 
negligible radiative losses, which makes the calculated energy loss directly 
comparable with the simulated energy deposit. The Geant4 simulation was 
checked against a different simulation [20] where energy losses are implemented 
at the microscopic level according to excitation and ionization cross sections. 
In particular, the energy deposited by relativistic electrons in the integrating 
sphere's nitrogen volume was compared. A systematic uncertainty of 2% on the 
simulation of the energy deposit was conservatively estimated, which includes 
the accuracy in the reference data of the Bethe-Bloch stopping power. This 
estimate is consistent with a detailed study of the simulation of the energy 
deposit presented in [30] . 

The emission of Cherenkov photons is a standard electromagnetic process 
in the Geant4 simulation, with an absolute yield determined by the index 
of refraction of the gas. Parameterizations of the index of refraction as a 
function of wavelength for nitrogen [31] [32], air [33] and helium [31] were 
implemented in the Geant4 simulation. Several checks were performed on the 
accuracy of these parameterizations. The nitrogen index of refraction taken 
from [31] is based on measurements between 468 nm and 2060 nm, while that 
taken from [32] is based on measurements between 145 nm and 270 nm. Both 
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parameterizations fit data between 230 nm and 490 nm [34J extremely well, 
differing only 0.2% from the index of refraction measured at 334 nm. The 
corresponding Cherenkov yields, integrated over the bandwidth of the 337 nm 
filter, were found to differ by less than 0.1%. For the refraction index of air, the 
uncertainty of Ciddor's parameterization [33J between 230 nm and 1700 nm 
is a few 10 -4 . From these studies, the systematic uncertainty on the absolute 
Cherenkov yield is estimated to be less than 0.1%. 



3.3 Simulation results 



Monte Carlo simulations of the experiment relevant for the fluorescence mea- 
surement were performed. Protons of 120 GeV energy were generated at 20 m 
distance from the apparatus, and tracked through the air volume and a beam 
tagging scintillator of the Test Beam line before reaching the AIRFLY ap- 
paratus. The beam profile was simulated according to the measured spread. 
For each simulated event, the energy deposited by the proton in the AIRFLY 
counters was recorded. Fluorescence and Cherenkov photons emitted along 
the proton's path inside the chamber were tracked and diffused within the 
integrating sphere. Photons exiting from the side port of the sphere may pass 
through the 337 nm filter and reach the PMT, where a photoelectron was 
detected according to the PMT efficiency An example of simulation is shown 
in Fig. where the number of detected photoelectrons from fluorescence and 
Cherenkov photons produced by 10 7 protons in air is given as a function of 
wavelength. 

The rates of single photoelectrons, expressed in units of photoelectrons per 
beam particle (pbp), are given in Tabled] as S^[ c and 5^ for the Fluorescence 
and Cherenkov simulated configurations, respectively. Results for both the 
diffuser plugs and the mirror plugs are reported. Also, the laser calibration 
constant, (Section I2.3[) . in units of photoelectrons per 337 nm photon 

(pp), was derived from a Monte Carlo simulation of the laser calibration. The 
statistical uncertainty on all simulation results is less than 1%. 





Sff (10- 4 pbp) 


s&° (io- 4 pbp) 


« r (io- 4 PP) 


diffuser plugs 


2.442 


10.53 


11.58 


mirror plugs 


2.131 


8.267 





Table 1 

Summary of Monte Carlo simulation results. 
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Fig. 5. Simulation of the AIRFLY experiment: number of detected photoelectrons 
from fluorescence (full line) and Cherenkov (dot-dashed line) photons produced by 
10 7 protons in air. 



4 Data analysis 



The rate of single photoelectrons at the PMT in time coincidence with a 
single beam particle was measured to determine the fluorescence yield. The 
selection of single beam particles was strongly based on the properties of 
the Cherenkov counter tagging the beam at the exit of the chamber. For 
illustration, a spectrum of the Cherenkov counter, obtained by integrating 
the charge around the known positions of the bunches (Fig. [2]), is shown in 
Fig. [6] for a data run where the beam was purposely misaligned. The peak (3) 
around 2300 ADC counts, corresponding to a single beam particle passing 
through the rod, is well separated from the pedestal (1) and from the signal 
for two beam particles (4). Also, a peak (2) is present corresponding to beam 
particles hitting the photomultiplier tube's window outside the 1 cm diameter 
acceptance of the rod. Thus, the Cherenkov counter provided a very clean 
single particle tagging, in both charge resolution and space. 

A single particle candidate in the Cherenkov counter was validated by a time 
coincident detection in the scintillator disk counter at the entrance of the 
AIRFLY chamber. In addition, each candidate in a train was required to be 
well isolated in a window of ±40 ns, with no other single particle signal in 
the Cherenkov and scintillator disk counters. In the same time window, veto 
counters were required to have less than 10% of their single particle signal in 
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Fig. 6. Spectrum of the Cherenkov counter: (1) pedestal, truncated for better vi- 
sualization; (2) signal corresponding to beam particles hitting the photomultiplier 
tube's window outside the 1 cm diameter acceptance of the rod; (3) signal corre- 
sponding to a single beam particle passing through the rod; (4) signal corresponding 
to two beam particles. 

order to reject off-axis beam particles. Lastly, the whole train was rejected if 
more than four single particle candidates were found in the Cherenkov counter 
and/or in the scintillator disk counter. Based on these strict criteria, about 
25% of all trains were selected, with a multiplicity of about 1.1 beam particle 
per train. To illustrate the quality of the single particle selection, the correla- 
tion between the signals of the Cherenkov and the scintillator disk counters is 
shown in Fig. [7] for one run. Events falling in the square are selected as single 
beam particles. Events accumulating in the upper right part of the figure, well 
separated from the selection area, correspond to two beam particles passing in 
the counters, and amount to only w 1% of all events. Events in the band with 
scintillator disk signal less than « 1000 ADC counts, consistent with a single 
beam particle in the Cherenkov counter, are produced by particles hitting the 
less efficient border of the scintillator disk, and amount to only « 0.5% of all 
events. From these data, the contamination in the selected sample from events 
of multiplicity higher than one is estimated to be negligible. 

Then, a single photoelectron in time coincidence with the beam particle was 
searched for in the PMT trace. The distribution of the ADC charge integral 
corresponding to peaks in the PMT trace is shown in Fig. [8] for one data run 
taken in the Cherenkov configuration. To count a single photoelectron, the 
charge integral was required to be between 300 and 3000 ADC counts. 
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Fig. 7. Correlation between the signals of the Cherenkov counter and the scintillator 
disk counter. Events within the square were selected as a single beam particles. 
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Fig. 8. The distribution of the ADC charge integral corresponding to peaks in the 
PMT trace. A single photoelectron signal is evident. 
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Fig. 9. Distribution of the arrival times of single photoelectrons relative to the arrival 
time of the beam particle. 

The arrival time of the selected photoelectrons relative to the beam particle, 
At, is shown in Fig. |9j A clear signal is observed around At = extending up 
to 30 ns. The exponential tail in the time distribution is due to the integrating 
sphere, since photons experience multiple diffusions before detection. Notice 
that a small background of single photoelectrons from accidentals is present 
before and after the main peak. The total number of single photoelectrons in 
a data taking run was calculated by integrating the time distribution in the 
interval [-4 ns, 32 ns], and subtracting the accidentals background as estimated 
from counts before -4 ns. 

The rate of single photoelectrons, in units of pbp, was then obtained by di- 
viding the total number of single photoelectrons by the total number of single 
beam particles selected in the run. Results from several runs taken in the same 
configuration were averaged to improve the statistical uncertainty. In the fol- 
lowing, the rates measured in a given gas will be denoted by i?ff and Rq^ for 
the Fluorescence and Cherenkov configurations, respectively. 

Results are summarized in Tables |2] to HJ where the quoted uncertainties are 
statistical only. The rates of single photoelectrons measured with the diffuser 
plugs are reported in Table [2] Also, measurements with the chamber under 
vacuum yielded: 

R£ c = (0-48 ± 0.09) ■ 10~ 4 pbp, (2) 
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= (3.11 ±0.21) • 1(T 4 pbp. 



(3) 



The rates of single photoelectrons measured with the mirror plugs are reported 
in Table El A set of runs was dedicated to fluorescence measurements with 



pN 2 
K Fl 

(10- 4 pbp) 


pair 

n m 
(10- 4 pbp) 


pHe 
K F\ 

(10- 4 pbp) 


pN 2 
^Ch 

(10- 4 pbp) 


pair 
^Ch 

(10- 4 pbp) 


pHe 
^Ch 

(10~ 4 pbp) 


20.05 ±0.11 


3.22 ± 0.07 


1.11 ±0.06 


32.89 ±0.15 


16.07 ±0.21 


3.78 ±0.11 


Table 2 

Summary of the rates of single photoelectrons measured with the diffuser plugs. 




(1(T 4 pbp) 


pHe 

(10- 4 pbp) 


(10- 4 pbp) 


pHe 
^Ch 

(10- 4 pbp) 






17.66 ±0.19 


1.19 ±0.09 


25.78 ±0.21 


1.05 ± 0.09 



Table 3 



Summary of the rates of single photoelectrons measured with the mirror plugs. 

an associated UV laser calibration (Section 12. 3f) . The diffuser plugs were used 
for these measurements. Results are summarized in Table HI For the measure- 
ment of the laser calibration constant, S'laser, single photoelectrons, selected 
with the same criteria as the fluorescence measurement, were counted in time 
coincidence with the laser shots. 



pN 2 

(10- 4 pbp) 


pHc 
K F1 

(10- 4 pbp) 


"S'laser 

(10- 4 pp) 


19.31 ±0.14 


1.20 ±0.08 


10.70 ±0.11 



Table 4 

Summary of measured fluorescence rates with an associated UV laser calibration. 
.S'laser is the measured laser calibration constant. 



5 Fluorescence yield measurement method 

The first step in the measurement of the fluorescence yield is the determination 
of the fluorescence rate in the 337 nm band, S Fl s , for the atmospheric gas under 
study. The rate, R^ s , measured in the Fluorescence configuration is given by: 

Tjgas _ ngas . r>gas / ,\ 

K Fl - b Fl + B Fl , (4) 

where Bp" 5 accounts for the background. Sources of background could include 
particles - from the beam halo or from secondary interactions - producing 
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Cherenkov and scintillation light in the PMT glass, in the interference filter 
and in the PTFE diffuser, or producing fluorescence in the small air volume 
between the PMT and the filter. The background can be subtracted with 
high accuracy by combining measurements in nitrogen and air. In fact, since 
the production of secondary particles - mainly electrons - is practically the 
same in nitrogen and air, and particle interactions in the apparatus and its 
surroundings do not depend on the gas filling the chamber, Bpf = Bp{ = B-p\ 
holds. Thus, backgrounds cancel in the difference, Ai?pi, between the rates 
measured in nitrogen and air: 



with tn 2 given by Eq. [U The fluorescence rate in nitrogen, S* F1 2 , or equivalently 
the rate in air, Spf = Spf /r^ 2 , is obtained from Eq. |5]and the measured Ai?pi- 

The measured fluorescence rate is converted into a fluorescence yield by two 
independent calibration methods. The first method is based on the calibrated 
337 nm light of a nitrogen laser. The procedure to derive the laser calibration 
constant, Si^, is described in detail in Section 12.31 The second method is 
based on the known amount of Cherenkov light produced by the beam par- 
ticle. Both Cherenkov and fluorescence emission contribute in the Cherenkov 
configuration, and the measured rate, Rq£ is given by: 



where is the rate due to Cherenkov light and Bq£ accounts for background 
originating from interactions of the beam particles in the exit port plug. To 
separate the Cherenkov contribution, the difference, AR^_ Fl , between the 
rates measured in the same gas with the Cherenkov and the Fluorescence 
configuration (Eqs. H] and ED is taken: 

A JD9 as _ vgas rjgas _ qgas . pgas /y\ 
^•n-Ch-Fl — -"-Ch n Fl ~ °Ch n Ch ■ \ 1 ) 

When diffuser plugs are used, Bq^ is dominated by Cherenkov light produced 
in the exit port plug and diffused back into the integrating sphere. When 
mirror plugs are used, a much smaller background from transition radiation 
by the thin aluminized mylar foil is expected. In both cases, the background 
does not depend on the particular gas filling the chamber. Secondary parti- 
cles, produced by the interaction of the beam with the exit port plug, are 
another potential source of background, since they could induce emission of 
fluorescence light in the gas and of Cherenkov light in the sphere's diffuser. 
This background contribution, which would depend on the gas, was found 
to be negligible from both Monte Carlo simulations and cross-checks with 



AR Fl 




(5) 




(6) 
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the data (Section l6.3p . Thus, Bq£ = Bch holds, and the background can be 
effectively eliminated by an appropriate combination of measurements with 
different gases. In particular, measurements with helium, which accompanied 
every measurement with atmospheric gases, are used to measure the difference: 

A f? N 2 A F? He — Q N 2 QHe /n\ 

^-"■Ch-Fl AA - ft Ch-Fl — °Ch °Ch- \°) 



The Cherenkov rate in nitrogen is derived from Eq. [HJ accounting for the small 
Cherenkov signal in helium (Sq^/Sq^ = 1%). 

The fluorescence rate, Sp{, and the calibration rate, S ca \ = Sq^ or S\ aseV} are 
then combined in the ratio Sp{ / S ca \. Notice that the absolute photon detection 
efficiency of the experiment cancels in the ratio. Thus, systematic uncertainties 
of the fluorescence yield associated with the efficiency of the integrating sphere, 
the transmission of the interference filter, and most notably the PMT quantum 
efficiency are significantly reduced. An equivalent ratio of S^ c to S^p = 
or Sj^ger is obtained from the Monte Carlo simulation of the experiment 
(Section 13.31) using a nominal fluorescence yield = 5.15 photons/Me V. 

The fluorescence yield of the 337 nm band in air is then measured from: 



( Sn 2 /r N2 )/ S ca u 

CMC / Of 
°F1 / °cal 



veal 

I 337 ~ oMC/QMC J 337 • \ J ) 



6 Fluorescence yield with Cherenkov calibration 



In this Section, measurements of the fluorescence yield based on Cherenkov 
light for absolute calibration are presented. Two measurements were per- 
formed, with a slightly different apparatus (Section [2,2]) . In the first, Cherenkov 
light emitted by the beam particle was diffused back into the integrating sphere 
by the diffuser mounted on the exit port plug. In the second measurement, 
the Cherenkov light was reflected back into the sphere by the aluminized my- 
lar mirror stretched over the exit port plug. In the following, details of the 
measurements are given. 



6. 1 Measurement with the diffuser plugs 



The analysis procedure detailed in Section [5] was used to measure the fluores- 
cence yield, using data collected with the diffuser plugs mounted on the top 
and exit ports of the integrating sphere (Table |2]). 



21 



The fluorescence rate in nitrogen is derived from Eq. [5] and the measured 
difference: 



AR m = Rpf - R$ = (16.83 ± 0.13) • 10" 4 pbp, (10) 
giving: 

5N2 = ( lg 44 ± 0.15) . 10- 4 pbp. (11) 

Also, the fluorescence background is estimated from Eq. [4] as: 

B Fl = (0.61 ± 0.07) • 10~ 4 pbp, (12) 



which is only 3% of the measured rate R 



N 2 
Fl • 



The Cherenkov rate is derived from the difference (Eq. [8]): 

A^_ F1 - A^e_ pi = s *h _ s nc = (1Q 1? ± 23) . 1Q -4 pbp (13) 

Notice that Rpf = (1.11 ± 0.06) • 10~ 4 pbp is larger than the fluorescence 
background of Eq. [121 indicating a nonzero Spf. Helium measurements were 
found to be very stable over the beam data taking, period, suggesting a residual 
nitrogen contamination in the helium gas bottlqjj. This small fluorescence-like 
signal in helium does not affect Eq. [T3J since S^f cancels in Ai?^_ F1 . After 
correcting Eq. [13] for the small Cherenkov signal in helium, the Cherenkov 
rate in nitrogen is obtained: 

S£g = (10.27 ± 0.23) • 10~ 4 pbp. (14) 

The Cherenkov background is estimated as: 

B Ch = AR%l_ m - SJ& = (2.57 ± 0.13) • 10~ 4 pbp. (15) 



An independent measurement of the Cherenkov signal was obtained from runs 
performed in the Cherenkov configuration with air filling the chamber. Follow- 
ing the same procedure used to derive the signal in nitrogen, the Cherenkov 
rate in air was measured to be: 

= (10.28 ± 0.26) ■ 10" 4 pbp. (16) 



3 Nitrogen impurities in helium are known to produce sizable fluorescence emis- 
sion [35] . Also, we observed a strong signal when introducing in our chamber a few 
hPa of nitrogen gas in helium at atmospheric pressure. 
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Due to the slightly different index of refraction, the Cherenkov yield in ni- 
trogen is 1.7% greater than the one in air, well within the uncertainty of the 
measurements. In the Monte Carlo simulation, the pressure and temperature 
of the air volume were chosen to match the average index of refraction of the 
air and nitrogen measurements. Thus, the best estimate of the Cherenkov rate 
is given by the average of Eq. [ID and Eq. [T6j 

= (10.27 ± 0.19) • KT 4 pbp, (17) 

in very good agreement with the Monte Carlo simulation (S^ = 10.53 • 
10 -4 pbp, Tabled]). This agreement, while not essential for the fluorescence 
yield measurement since the absolute photon detection efficiency cancels in 
the ratio of the fluorescence to Cherenkov rate, confirms the reliability of the 
simulation in reproducing the overall efficiency of the experiment. 

Equation O with S CSL \ = Sq£ and S^p = Sq^, is then used to derive the fluo- 
rescence yield. From the measured fluorescence and Cherenkov rates (Eqs. [TT1 
and [IT]), and the corresponding Monte Carlo rates of Table [I] the fluorescence 
yield of the 337 nm band in air is found to be: 

FgChDiff = 5 64 ± o.i2 photons/MeV. (18) 
6.2 Measurement with the mirror plugs 

Data collected with the mirror plugs mounted on the top and exit ports of 
the integrating sphere (Table [3]) were also analyzed following the procedure 
described in Section [5j 

The fluorescence rate in nitrogen is obtained by subtracting the background 
of Eq. [T^l from the measured signal in the Fluorescence configuration: 

Sgi = #N 2 _ = ( 1? Q5 ± 21 ) . 1Q -4 pbp ( 1Q ) 

The Cherenkov rate is derived from the difference (Eq. [8]): 

Ai?^ _ Ai?ge_ Fi = Sgg - S% = (8.26 ± 0.31) ■ 10~ 4 pbp. (20) 

After correcting Eq.[20]for the small Cherenkov signal in helium, the Cherenkov 
rate in nitrogen is obtained: 

= (8-34 ± 0.31) • HT 4 pbp, (21) 
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in very good agreement with the Monte Carlo simulation (S^f = 8.27 ■ 
10~ 4 pbp, Table [1]). This agreement gives further confidence on the quality 
of the simulation. 

The fluorescence yield is then derived from Eq. with 5 ca j = S c ^ and S^p = 
S^h- From the measured fluorescence and Cherenkov rates (Eqs. [T9l and l2Tj) , 
and the corresponding Monte Carlo rates of Table [U the fluorescence yield of 
the 337 nm band in air is found to be: 

yChMirr = g^g ± q 35 photons / MeV> (22) 



in very good agreement with Y^ 01 ^ (Eq. [T8|) . 

Notice that the absolute photon detection efficiency is ~ 20% smaller with 
the mirror plugs (Table [1]). Also, there is a significant difference between 
the Cherenkov background with the mirror plugs, estimated as Ai?^£_ F1 = 
(-0.14±0.13)-10" 4 pbp, and that with the diffuser plugs, B Ch = 2.57-10" 4 pbp 
(Eq. [To^l . The agreement between y 3 gt; Mirr and Y^ DlS confirms the reliability 
of the Cherenkov calibration method. 



6.3 Cross-checks and systematic uncertainties 



Systematic uncertainties are largely correlated in the two Cherenkov calibra- 
tion measurements, which use the same apparatus - except for the top and 
exit port plugs of the integrating sphere - and data analysis technique. In the 
following, we will refer to Y^ when cross-checks and systematic uncertainties 
apply to both K^ Diff and Yg^ MiTI . 

Fluorescence of the diffuser material could potentially bias the fluorescence 
yield measurement. Negligible fluorescence induced by UV light, advertised 
by the GORE™ DRP® manufacturer, was confirmed by our measurements 
with a deuterium lamp. In addition, fluorescence produced by charged parti- 
cles traversing the diffuser was also measured with an 18 /zCi 241 Am source 
placed at the entrance port of the sphere. The chamber was evacuated, letting 
alpha particles hit the inner surface of the sphere. In these conditions, the 
PMT single photoelectron rate, ~ 800 Hz, was found to be independent of the 
presence of the 241 Am source. The PMT rate increased to 20 kHz when nitro- 
gen at atmospheric pressure filled the chamber, due to fluorescence emission 
from the nitrogen molecules excited by the alpha particles. From these results, 
fluorescence of the diffuser induced by charged particles, if any, can be safely 
considered negligible as far as the Y337 measurement is concerned. 

In the Fluorescence configuration, Cherenkov light produced by the beam 
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particles is dumped on the black UV absorbing material covering the inner 
surface of the chamber. A fraction of this light could get back into the sphere 
through its open exit port, introducing a bias on Y337. To investigate this 
potential systematic effect, a 337 nm laser beam was shot along the beam 
line into the Cherenkov dump, and the corresponding light intensity at the 
PMT port was measured. The light intensity was found to be < 10 -3 of that 
measured when the laser hit the sphere's diffusive wall. Since the nitrogen 
fluorescence light collected by the sphere is about twice the Cherenkov light 
(Sections 16.11 and 16.21) . the effect of Cherenkov light bouncing back into the 
sphere from the Cherenkov light dump is negligible. 

Several checks of the data analysis and background subtraction method were 
performed. K537 was found not to depend on the value of the threshold chosen 
to select single photoelectrons. This result was expected, since the same thresh- 
old is used in both fluorescence and Cherenkov measurements. A potential bias 
due to the slightly different time evolution of the Cherenkov and fluorescence 
signals was studied by changing the time window for photon counting between 
24 and 40 ns, which produced differences in Y337 of less than 1%. Much stricter 
selection criteria for isolated particles, which reduced the selected sample by a 
factor of two, resulted in a 1% shift of K^. An independently developed data 
selection analysis, with a different peak finding algorithm and background re- 
jection criteria, was also performed, which changed Y^j by 0.7%. Potential 
biases in the determination of the fluorescence and Cherenkov backgrounds 
were also investigated. For Y^ DlS , measurements with the chamber evacuated 
provide an independent estimate of B F y and B C h- Rp\ c = (0.48±0.09)-10~ 4 pbp 
(Eq. [2]) is consistent within uncertainties with the fluorescence background of 
Eq. [121 confirming the Monte Carlo expectation that background from secon- 
daries produced in the chamber gas is negligible. An estimate of Bch is derived 
from AR}%_ Fl = (2.63 ± 0.23) • 10~ 4 pbp, which is consistent with Eq. QU An- 
other check of B Ch is provided by i^-i2g£ = (16.82 ±0.26) ■ 10" 4 pbp, where 
the Cherenkov background in nitrogen and air should cancel in the difference. 
The result is consistent with Eq. [10] within their uncertainties, indicating that 
Cherenkov background does not depend on the chamber gas, and fluorescence 
produced by secondaries from the exit port plug is negligible, as expected 
from Monte Carlo simulations. For , the same fluorescence background 

(Eq. [T2|) of the diffuser plugs measurement is used. This choice is justified by 
the stability of the helium measurements (Tables [2] and [3]), performed routinely 
to monitor the background level. If the observed difference (0.08 ± 0.11 pbp) 
between the two i?pf rates is due to a change of the fluorescence background 
level, the effect on Y^j DlS would only be of 0.5%. All performed checks resulted 
in shifts of Y^fy smaller than its statistical uncertainty, and a 1% systematic 
uncertainty is conservatively assigned to the data selection and background 
subtraction methods. 

Recent studies [35] indicate that the absolute value of r Na may depend on the 
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presence of very small amounts of contaminants in the nitrogen gas. Thus, a 
bias in the fluorescence yield could be introduced by changes in the nitrogen 
gas composition during the measurements. To study the stability of r Na , we 
performed measurements with different gas bottles, with different conditions 
of gas flow in the chamber and for extended periods of time. The observed vari- 
ations are included in the 1% uncertainty quoted in Eq. [IJ Also, the rates R F f 
measured during several weeks at the Test Beam did not show any noticeable 
change beyond statistical fluctuations. From these studies, a 1% systematic 
uncertainty is assigned to for the nitrogen to air fluorescence ratio. 

The absolute value of the photon detection efficiency, including the peak value 
of the integrating sphere's efficiency, of the 337 nm filter transmission and of 
the PMT quantum and collection efficiency, cancels in the ratio of fluorescence 
to Cherenkov signal, and its contribution to the systematic uncertainty of 
is negligible. 

A difference in the sphere's efficiency between the Fluorescence and Cherenkov 
configurations may introduce a bias. This was studied by shooting a 337 nm 
laser beam on different points of the sphere's inner surface, and switching 
between the Fluorescence and Cherenkov configurations. The maximum shift 
in the light intensity measured at the side port was 0.9%, which is taken as 
an estimate of the systematic uncertainty on Y^j associated to the sphere's 
efficiency. 

The uncertainty associated to the wavelength dependence of the detection 
efficiency must also be evaluated, due to the wide wavelength range of the 
Cherenkov emission. In the Monte Carlo simulation of the experiment, the 
wavelength dependence of the diffuser is included in each diffusion bounce 
inside the sphere (Section 13. 1 f) . We also measured the overall efficiency of the 
integrating sphere between 300 and 400 nm by detecting at the side port the 
light spectrum of a deuterium lamp with an Oriel MS257™ spectrograph. 
This efficiency was then used in a simulation to weight photons exiting the 
side port. The two simulation approaches yielded a 1% shift of Y£$, which 
is taken as an estimate of the systematic uncertainty. The wavelength depen- 
dence of the PMT quantum efficiency was varied according to its uncertainty, 
and a corresponding 1% systematic uncertainty is assigned. The wavelength 
and angular dependence transmission of the 337 nm interference filter was 
measured by two independent groups; when used in the Monte Carlo simula- 
tion, the difference between the two measurements resulted in a 2% shift of 
Kj37, which is taken as systematic uncertainty. The uncertainty in the relative 
intensities of the fluorescence bands [T2] has a negligible effect, since the signal 
is dominated by the 337 nm band (Fig. [5]). 

From the studies of Section |3~2| a 2% systematic uncertainty is assigned to the 
simulation of the energy deposit, while the uncertainty of the Cherenkov yield 
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from the gas index of refraction is considered negligible. The limited statistics 
of the Monte Carlo simulations contribute 1% to the uncertainty. 

The different contributions to the systematic uncertainty are summarized in 
Table EJ Taking these contributions as uncorrelated, a total systematic un- 
certainty of 3.7% on is obtained. Notice that all contributions should be 
considered fully correlated between Y^ DlS and Y^ Mm , with the exception of 
the Monte Carlo statistics since separate simulations were performed for the 
two setups. 



data selection and background subtraction 


1.0% 




1.0% 


integrating sphere efficiency 


0.9% 


integrating sphere wavelength dependence 


1.0% 


PMT quantum efficiency 


1.0% 


filter transmittance 


2.0% 


simulation of energy deposit 


2.0% 


Monte Carlo statistics 


1.0% 


Total 


3.7% 



Table 5 

Systematic uncertainties on the fluorescence yield measured with Cherenkov cali- 
bration, y 3 ^ Diff and Yg^ Mirr . 
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7 Fluorescence yield with laser calibration 



The fluorescence yield measurement reported in this Section is largely inde- 
pendent of those based on Cherenkov calibration, since a nitrogen laser (Sec- 
tion |23]) was used for the absolute calibration. Data presented in Table H] were 
analyzed with the procedure of Sectional 

The fluorescence rate in nitrogen is obtained by subtracting the background 
of Eq. [12] from the rate measured in the Fluorescence configuration: 

£g> = R^f - B Fl = (18.70 ± 0.16) ■ 10~ 4 pbp. (23) 



For the absolute calibration, the laser calibration constant (Table H]) 

S laser = (10.70 ± 0.11) • 10" 4 pp, (24) 



was obtained from several calibration runs performed routinely at the end of 
each day of data taking. 

The fluorescence yield is then derived from Eq. [H] with S ca _i = Siaser and 
ScaF = ^laser- From the measured fluorescence rate (Eq. [23]), the laser calibra- 
tion constant (Eq. [241 . and the corresponding Monte Carlo rates of Table [1] 
the fluorescence yield of the 337 nm band in air is found to be: 

Y^ 7 eT = 5.73 ± 0.08 photons/MeV, (25) 



in very good agreement with the Cherenkov calibrated measurements Y^ DlS 
(Eq. USD and Kg!? Mirr (Eq. [22]). 

Notice that the fluorescence rate (Eq. 1231) is somewhat smaller than the rate 
measured in the Cherenkov calibration measurement (Eq. [TT1) . In fact, after 
completing the measurement with laser calibration, the apparatus was dis- 
assembled, the integrating sphere was cleaned and new diffuser plugs were 
installed, which resulted in a slightly higher efficiency of the sphere in the 
subsequent measurement with Cherenkov calibration. Since the efficiency of 
the sphere cancels out in Eq. [9] the fluorescence measurement is not affected 
and the agreement between Y^j er and the measurements with Cherenkov cal- 
ibration confirms the reliability of the calibration methods. 
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7.1 Cross-checks and systematic uncertainties 



Several of the systematic uncertainties evaluated in Section 16.31 are common 
to Y^ er , since the same apparatus and data analysis techniques were used. In 
particular, systematic uncertainties from the data selection and background 
subtraction, the nitrogen to air fluorescence ratio, the integrating sphere's 
efficiency, and the simulation of the energy deposit are fully correlated between 
^337 er an d ^337- On the other hand, systematic uncertainties associated to 
the wavelength dependence of the integrating sphere's efficiency, to the PMT 
quantum efficiency, and to the interference filter were found to be negligible, 
since the nitrogen laser calibration light, unlike Cherenkov light, reproduces 
exactly the 337 nm fluorescence band. 

Checks on the background subtraction specific to the laser calibration mea- 
surement included the stability of the helium measurements (Tables [2] and H]), 
performed routinely to monitor the background level. If the observed differ- 
ence (0.09 ± 0.10 pbp) between the two Rpf rates is attributed to a change of 
the fluorescence background level, the effect on Y^, eT would only be of 0.5%. 
The background contribution from accidentals in the determination of Si^er 
was negligible, and does not contribute to the systematic uncertainty. 

Systematic uncertainties specific to the laser calibration procedure included 
the knowledge of the geometry of the chamber. Since the fluorescence signal is 
proportional to the path length of the beam particle in the integrating sphere, 
Sp\ c depends on the geometry introduced in the Monte Carlo simulation. 
Due to the mechanical assembly of the integrating sphere, which includes in 
its interior the soft PTFE layer of 0.5 mm thickness, the uncertainty on the 
internal size of the sphere is conservatively taken as 0.5 mm, which corresponds 
to a 0.3% systematic uncertainty on Kj'ff*. Notice that this uncertainty is 
absent in the Cherenkov calibration measurements of Section [6], since both 
fluorescence and Cherenkov light are proportional to the path length of the 
beam particle, which cancels in their ratio. 

The absolute calibration of the laser probe and the efficiency of the laser cali- 
bration sphere introduce systematic uncertainties specifically associated with 
the laser calibration procedure. For the absolute calibration of the RjP-465 
probe, the 5% uncertainty quoted by the manufacturer was taken. The repro- 
ducibility of the absolute calibration was checked by comparing two RjP-465 
probes, which were found to agree within 2.6%, well within the quoted system- 
atic uncertainty. For the efficiency of the laser calibration sphere (Section 12.3ft . 
the 0.5% uncertainty associated to the linearity of the probe, combined with 
a 0.6% uncertainty estimated from the reproducibility of the measurement, 
gave a 0.8% systematic uncertainty on Y^ er . 
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Other checks include the stability of S^ser and measurements with different 
particle beams. The laser calibration constant was measured shooting the laser 
beam on different spots in the integrating sphere, and repeating the measure- 
ment after disassembling and remounting the laser system. In all of these 
checks, S'laser changed by less than 1%. Measurements of the fluorescence yield 
with beams of 32 GeV pions and 8 GeV positrons were also performed during 
the test beam period dedicated to the laser calibration. The multiplicity of 
these secondary beams was significantly lower than the primary proton beam, 
and the corresponding fluorescence yield measurements are not statistically 
competitive. Nevertheless, they provide a useful cross-check of the indepen- 
dence of the fluorescence yield to the type of relativistic particle originating the 
energy deposit. Indeed, measurements performed with these secondary beams 
were consistent within their 6% statistical uncertainty with the fluorescence 
yield of Eq. E3 

The different contributions to the systematic uncertainty are summarized in 
Table [61 Taking these contributions as uncorrelated, a total systematic uncer- 
tainty of 5.8% on Kjff er i s obtained. 



data selection and background subtraction 


1.0% 




1.0% 


integrating sphere efficiency 


0.9% 


geometry 


0.3% 


laser probe calibration 


5.0% 


calibration sphere transmission 


0.8% 


simulation of energy deposit 


2.0% 


Monte Carlo statistics 


1.0% 


Total 


5.8% 



Table 6 

Systematic uncertainties on the fluorescence yield measured with laser calibration, 

\^laser 
I 337 ■ 



8 Combined fluorescence yield measurement 



The fluorescence yield measurements presented in Sections [6] and [7] - F 3 37 Dlff , 
y 3 37 Mirr and y^f| er - are found to be in good agreement within their uncertain- 
ties. They are based on very different calibration light sources - photons from 
Cherenkov emission or from a nitrogen laser - whose corresponding systematic 
uncertainties are largely uncorrelated. Thus, an appropriate weighted average 
of these measurements yields an improved result. 
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First, the measured fluorescence yields were extrapolated to a nominal pres- 
sure of 1013 hPa and temperature of 293 K. This procedure, which used our 
measured p' air (337) = 15.89 hPa [12] and a 337 = —0.36 [32], shifted the mea- 
surements by less than 0.5%. Then, a weighted average was performed, tak- 
ing into account the correlation between uncertainties. In particular, system- 
atic uncertainties associated to the integrating sphere wavelength dependence, 
the PMT quantum efficiency, the filter transmittance, the geometry, the laser 
probe calibration, and the calibration sphere transmission were taken as un- 
corrected between the Cherenkov and the laser calibration measurements. 
From the weighted average, the fluorescence yield of the 337 nm band in air 
at 1013 hPa and 293 K was found to be: 

r 337 = 5.61 ± 0.0Q stat ± 0.21 syst photons/MeV, (26) 

where the statistical and systematic uncertainties are quoted separately. 

Of the recent fluorescence yield experiments, only a few [6] [10] have measured 
the yield of the 337 nm band in air and can be directly compared with the 
AIRFLY measurement of Eq. [2J3 Most of the experiments [5] [7] |8J [9J report 
a fluorescence yield integrated over the spectrum between « 300 to 400 nm. 
To compare with these experiments, we convert the integrated yield into a 
yield of the 337 nm band by using the spectrum measured by AIRFLY [12]. 
Results for Y337 are compared in Fig. [IOJ The AIRFLY measurement is com- 
patible with previous measurements and presents the smallest uncertainty. It 
is also compatible with the detailed study by [27], where it is argued that 
some of these measurements should be corrected for a systematic bias in their 
calculation of the energy deposit. 



9 Conclusions 

We have performed a precise measurement of the absolute fluorescence yield 
of the 337 nm band in air relevant for UHECR experiments. The fluorescence 
emission was calibrated by a known light source - Cherenkov emission from 
the beam particle or a calibrated nitrogen laser - measured in the same ap- 
paratus. With this novel experimental method, the impact of the photomulti- 
plier detection efficiency was minimized, reducing significantly the systematic 
uncertainty. Two independent calibrations provided consistent results, and a 
total uncertainty of 4% on the absolute fluorescence yield was achieved. 

The AIRFLY measurements - the absolute yield reported here and the pres- 
sure, temperature and humidity dependence of the fluorescence spectrum [12] 
[T5] - provide the most comprehensive and precise parameterization of the 
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Fig. 10. Experimental results on Y337. For some experiments, the fluorescence yield 
of the 337 nm band is derived from the integrated yield measured between ~ 300 
to 400 nm (see text for details). 

fluorescence yield currently available. These measurements have direct im- 
plications for UHECR experiments which employ Fluorescence Detectors to 
determine the cosmic ray energy. For example, the absolute fluorescence yield 
of the 337 nm band reported here is 11% and 30% larger than that currently 
adopted by the Pierre Auger Observatory [T7] [37] and by the Telescope Ar- 
ray [18] [38], respectively. While the actual effect on the UHECR energy spec- 
trum also depends on the specific fluorescence spectrum adopted by these ex- 
periments, a downward shift of the energy scale by at least ~ 10% is implied 
by the AIRFLY result. At the same time, the uncertainty on the energy scale 
associated to the fluorescence yield, currently a major contribution [T7] [Tg] . 
will be reduced by a factor of about three. 

In principle, the experimental methods developed by AIRFLY could be fur- 
ther refined to improve the precision of the fluorescence yield. In particular, 
the 5% systematic uncertainty of the laser energy probe - the main systematic 
of the pulsed laser calibration method - may be reduced, or a continuous laser 
absolutely calibrated to 1-2% could be employed. However, the uncertainty on 
the energy scale of UHECR experiments is likely to be dominated by other 
contributions, including the absolute calibration of the fluorescence telescopes 
and the knowledge of the atmosphere. Thus, we expect the absolute fluores- 
cence yield measured by AIRFLY to remain a reference for the current and 
next generation of UHECR experiments. 
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